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SUMMARY 


A single-cell Eulerian photochemical air-quality simulation model was 
developed and validated for selected days of the 1976 St. Louis Regional Air 
Pollution Study (RAPS) data sets; parameter izations of variables in the model 
and validation studies using the model are discussed. Good agreement was 
obtained between measured and modeled concentrations of NO, CO, and NO 2 for all 
days simulated. The maximum concentration of O 3 was also predicted well. Pre- 
dicted species concentrations were relatively insensitive to small variations 
in CO and NO^ emissions and to the concentrations of species which are entrained 
as the mixed layer rises. Because of the rather small horizontal dimensions of 
the single cell (20 km by 20 km) , predicted species concentrations were sensi- 
tive to the advection of upwind ozone and its precursors. Chemical reactions 
were dominated by the chemistry of reactive hydrocarbons. Significant changes 
in the predicted concentrations of O 3 were also obtained when photolytic reac- 
tion rates were scaled from theoretical to atmospheric values to reflect the 
attenuation of solar radiation by particulates and aerosols. The uses of addi- 
tional types of data in the model formulation and validation procedures are 
discussed. 


INTRODUCTION 

The NASA Langley Research Center has programs to develop both in situ and 
remote sensors for environmental air-quality measurements on the regional to 
global scale. Initial testing of the sensors normally occurs over relatively 
small (local) spatial scales using several instrumented aircraft available at 
Wallops Flight Center, Jet Propulsion Laboratory, and Langley Research Center. 
From 1977 to 1979, a series of summer field-measurement programs, the South- 
eastern Virginia Urban Plume Studies (SEV-UPS) , were performed to evaluate 
several NASA remote sensors (refs. 1 and 2) . To assist in the interpreta- 
tion of these measurements, an urban-scale air-quality simulation model was 
developed. 

The model utilizes a single-volume element or cell centered over an urban 
area with a fixed horizontal length and a variable vertical height that depends 
on the height of the mixed layer. Transport through the cell and photochemistry 
within the cell are treated simultaneously. The present paper describes the 
validation of this model by using several data sets obtained during the 1 976 EPA 
Regional Air Pollution Study (RAPS) in St. Louis (ref. 3) . 

The general approach used in the single-cell air-quality model of Schere 
and Demerjian (ref. 4) was followed in the present study. The representation 
of the physical processes in the single-cell model is compatible with the type 
and quantity of data that are normally present in urban air-quality and emis- 
sions data bases. An Eulerian photochemical air-quality simulation model 
similar to the model developed by Schere and Demerjian was developed and tested 
by comparing the model results with measurements from selected days in the RAPS 
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data set. When the methods used to parameterize variables in the Schere and 
Demerjian model (ref. 4) were defined, those methods were also used to param- 
eterize this model. However, the parameter izations for several key variables 
were not described in sufficient detail for inclusion in this model. The vari- 
ables which were parameterized in this study, and consequently differ from 
Schere and Demerjian’s model (ref. 4), are as follows: chemical mechanism, 

emissions inventory, growth rate of the mixed layer, and background and entrain- 
ment concentrations. The chemical mechanism describes only gaseous chemical 
reactions. The model is applicable principally to clear-sky days; however, 
days with cloudy skies can be modeled if solar radiation measurements are avail- 
able. The justifications for using the present modified techniques and input 
parameters are discussed in the subsequent sections. As a result, the verifica- 
tion results of air-quality trends for St. Louis for July 23, 1976 (day 205), 
and August 3, 1976 (day 216), two days presented by Schere and Demerjian, 
are necessarily somewhat different from results presented previously by those 
authors. In addition, verification studies using additional clear-sky days 
in the RAPS data set, June 7, June 8, July 13, and August 13, 1976 (days 159, 
160, 195, and 226, respectively), are discussed. Day numbers instead of dates 
are used throughout the present paper. 


MATHEMATICAL MODELING APPROACH 

The basic working equations governing the prediction of species concentra- 
tions for the single cell model are as follows: 
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where 

^i,t concentration of species i at time t, ppmv-hr”*^ 

hourly averaged emissions flux of species i, ppmv-m-hr"^ 
height of mixing layer at time t, m 

^i,bg concentration of species i at time t-1 upwind of the modeling 

region, ppmv 

Ci,top concentration of species i above the mixing layer, ppmv 
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Ax horizontal length of cell, km 

U hourly averaged wind speed through cell, m-hr”^ 

Ri differential rate of change of species concentrations due to chemical 

reactions , p^xtiv-hr 

t time 

At time increment, hr 

The transport and chemistry of the atmospheric species are coupled in order to 
better simulate the air quality in an urban area as seen by the in situ sensors. 
The time rate of change of species concentrations as predicted by equation (1 ) 
is dependent upon each of the following terms on the right-hand side: emis- 

sions into the cell, the advection inflow through the cell, the dilution effect 
of the expanding volume of the cell as the mixed layer rises, the entrainment 
of any pollutant species aloft, and the production and/or destruction of 
species caused by chemical reactions within the volume. Dry deposition and sur- 
face scavenging are neglected because insufficient data exist to parameterize 
or validate these processes. Figure 1 shows the cell centered over St. Louis 
with a horizontal cell length x of 20 km. The Regional Air Monitoring (RAM) 
station locations used to obtain the RAPS data set (ref. 3) are indicated. 

The horizontal lengths of the cell are fixed in the beginning of the 
modeling period, the base of the cell is the ground, and the top of the cell is 
the mixing layer height, which normally increases during daylight hours. Tur- 
bulent mixing, which produces uniform species concentrations within the cell 
volume is assumed to occur; therefore, diffusion is neglected. The geographic 
position of the cell, which is arbitrarily defined, is constant throughout the 
simulation. The direction of advection through the cell is defined by the 
prevailing wind direction and, therefore, any change slightly during the simu- 
lation. Thus, horizontal fluxes through the model cell may not always enter the 
cell perpendicular to a cell face. 

Chemical species were separated into three groups based on the chemical 
lifetimes. The separation into groups was performed to reduce the order of the 
matrix that was integrated and thus reduce the time required for integration 
(ref. 5) . Concentrations of species in two groups were held constant while solu- 
tions for species concentrations in a third group were obtained. 

Group 1 contains species whose lifetimes of several days are much greater 
than the time step (1.5 minutes) for integration of species concentrations. 
Species in the group (which include CO 2 # H 2 O, O 2 , and inert gases) have low 
chemical reactivities. Concentrations of these species were held constant 
throughout the simulation. 

Group 2 contains species whose chemical lifetimes are several hours in 
length; diurnal variations in concentrations are routinely monitored. Species 
defined by the Environmental Protection Agency (ref. 6 ) as criteria pollutants 
(NO 2 , NO, O 3 , and CO*) and other stable compounds are members of the group. 
Chemical reactions couple species concentrations within this group to species 
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concentrations both within and outside of the group. Integration of the equa- 
tions describing the changes in concentrations 6 cj[/ 6 t is performed using the 
EPISODE algorithm (ref. 7). 

Group 3 contains species whose lifetimes are short, on the order of seconds 
or a few minutes. The species have low concentrations and high chemical reactiv- 
ities and include reaction intermediates and radicals such as OH, HO 2 , NO 3 , and 
RCX) 3 . The equations describing the concentration changes of species in the 
group reduce from differential to nonlinear algebraic form when the steady-state 
hypothesis ( 6 cj^/ 6 t = 0) is invoked. As a consequence of using the steady-state 
hypothesis, transport is neglected for species in this group. Since the alge- 
braic equations are coupled, they are solved iteratively until a stable solu- 
tion for all group 3 species concentrations is achieved. 

The sequence of computations used to obtain predicted species concentra- 
tions is as follows. Measured concentrations of the coupled species at 0530 
local standard time (LST) were used to calculate concentrations of the steady- 
state species (group 3) at 0530 LST. The complete set of species concentrations 
was then used to form differential expressions (eq. ( 1 )) for each of the coupled 
species. Concentrations of coupled species returned from EPISODE for the sub- 
sequent time step were then used as constants in the solution of steady-state 
species concentrations for the subsequent time step. The sequence of computa- 
tions was repeated for the remainder of the simulations. 


PHOTOCHEMICAL MECHANISM 

The Falls and Seinfeld mechanism (ref. 8 ) was used in this model to 
describe the chemical kinetics that occur in the urban troposphere. The 
mechanism describes the gaseous chemical reactions of ozone, nitrogen oxides, 
and hydrocarbons in an urban atmosphere. Heterogeneous chemistry is not 
treated. Rate constants for thermal reactions were adjusted hourly using 
Arrhenius kinetics to account for the diurnal temperature variations. 

Hydrocarbons in the mechanism are partitioned into classes based on bonding 
characteristics; that is, hydrocarbons are separated into alkanes, ethylene, 
other alkenes, formaldehyde, other aldehydes, organic nitrates, and alkoxy and 
alkyl peroxy radicals. Hydrocarbon chemistry generally dominates the prediction 
of diurnal concentrations of ozone in urban atmospheres. However, the RAPS data 
set does not include speciation of nonmethane hydrocarbons, which is a require- 
ment for the efficient utilization of the mechanism. The sensitivity of the 
mechanism to variations in nonmethane hydrocarbon concentrations is addressed in 
a subsequent section. 

Quantum yields for the photolysis of NO 2 in the laboratory were obtained 
from Jones and Bayes (ref. 9) and used as input in an algorithm developed by 
Schere and Demer jian (ref. 1 0) to canpute diurnal variations of all photolytic 
rate constants in the mechanism. The theoretical rate constants are computed 
for clear-sky conditions and are a function of actinic flux, date, time, and 
location of the modeled region. Photolytic rate constants in the atmosphere 
are reduced relative to their theoretical values because of the attenuation 
of solar radiation by aerosols and gases in the polluted atmosphere; Zafonte 
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et al. (ref. 11 ) parameterized the attenuation by introducing a zenith-angle 
dependence into the atmospheric value of the photolytic rate constant for ^ 2 * 
The rate constant for the photolysis of NO 2 used in this calculation kj ^02 

reduced relative to its theoretical value by using an empirical scaling factor 
parameterized for the St. Louis area by Schere and Demerjian (ref. 4) and based 
on the work of Zafonte et al. (ref. 11). The scaling factor is a function of 
the measured (UV to near-lR) total solar radiation (TSR ) , the solar zenith 
angle z, and the theoretical photolysis rate for NO 2 : 

JtN 02 (scaled) = [-0.0696 + 0.835(1 - oos z)](TSR) (64° ^ z < 90°) 


kj ,02 (scaled) = 0.4064(TSR) 


%02 (scaled) 


0.295 + 0.076 



(TSR) 


(47° S z < 64°) 
( 0 ° ^ z < 47°) 


Rate constants for other photolytic reactions in the mechanism were reduced by 
a scale factor of kj^Q^ (scaled) days chosen for simula- 
tions in this work were primarily clear-sky days; thus, smooth variations in the 
photolytic rate constants were obtained. 


SPECIES CONCENTRATIONS 

The rationale for selecting concentrations of species required to integrate 
equation ( 1 ) is discussed in this section. 

Concentrations of O 3 , CH 4 , NO, N 02 / and total hydrocarbons were measured 
at many stations within and surrounding the urban area shown in figure 1 
(ref. 3) . An hourly averaged concentration for each species was calculated 
by using measurements at all stations within the cell for a given hour. Average 
measured concentrations for 0530 LST, the starting time for simulations, were 
used to initialize the model. Average measured concentrations for later hours 
were compared with calculated concentrations. 

Methane and total hydrocarbon concentrations were measured in St. Louis, 
but since hydrogen speciation was not part of the archived RAPS data set, parti- 
tioning of nonmethane hydrocarbons into bonding classes present in the mechanism 
(ref. 8 ) was performed within the model. The partitioning process was based on 
the relative percentages of average ambient concentrations of hydrocarbons mea- 
sured in Los Angeles in 1965 (ref. 12). Table I shows the levels and percent- 
ages of hydrocarbons found in Los Angeles in 1965 along with scale factors and 
percentages of hydrocarbons used in the present study. The percentages of 
reactive hydrocarbons found in St. Louis were expected to be substantially less 
than those in Los Angeles. The scale factors used for all final simulations, 
which were varied in increments of 5 percent for each nonmethane hydrocarbon 
class, were determined by modeling days 205 and 216 simultaneously in order to 
obtain calculated ozone curves which agreed in both magnitude and shape with the 
measured curves. Meteorological conditions and ozone air-quality trends for 
days 205 and 216 were markedly different from each other; simultaneous attempts 
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to match the ozone for both days resulted in hydrocarbon scaling factors that 
were more general than those which would have been obtained if each day had been 
considered separately. The effect of varying the scale factors on the predic- 
tion of the magnitude of the ozone concentration for day 21 6 is shown in fig- 
ure 2. Note that the hours are denoted by 5^ Ir etc., instead of 0500, 0700, etc. 
The measured ozone concentration is denoted by circles; the shaded area repre- 
sents a range of one standard deviation about the mean for ozone concentrations 
measured at all stations within the model cell at each hour- A change in the 
scale factors had no effect on the computed CO curve and only minor effects on 
the computed NO and NO 2 curves. While the scale factors used in this study are 
not definitive for the St. Louis study, they provided the best agreement of 
modeled and measured results for the days simulated and for the scale factors 
considered. More recent measurements of hydrocarbon speciation in major 
U.S. cities (including St. Louis) show reactive hydrocarbon partitioning coeffi- 
cients which are quite similar to those used in this model (ref. 13). 

Nonzero concentrations for all species are required to initiate the 
chemical-reaction portion of the model. Initial concentrations of species not 
measured but included in the mechanism were estimated by using the values of 
Duewer (ref. 14) when available. When not measured or available from Duewer, 
initial species concentrations inside the cell were set equal to the concentra- 
tions upwind of the cell (See table II.) 

Wind-direction measurements were used to select stations that measured 
hourly concentrations of coupled species upwind of the cell An axis 

system containing 45^ arcs was centered over the modeling region. The midpoints 
of the arcs correspond to N(O^), NE(45^), E(90^), etc. Stations were selected 
when their locations corresponded to the measured wind direction. 

Measurements of upwind concentrations of coupled species were not available 
for all hours simulated. Since hourly nonzero values are required to describe 
the transport through the cell, a sensitivity study was devised to estimate 
upwind concentrations. Measured initial concentrations inside the cell, upwind 
concentrations, concentration estimates from Duewer (ref. 14), and first-guess 
concentrations for coupled species were used to simulate diurnal concentrations 
for days 205 and 216. The first-guess concentrations were then successively 
varied, one species at a time, and the effects of a single concentration varia- 
tion on the simulated concentrations of O 3 , NO, NO 2 , and CO for days 205 and 216 
were examined. Unknown species concentrations (such as aldehydes) were reduced 
to the point where small variations in concentrations had no effect on computed 
ozone concentrations. 

An example of the effects of varying first-guess concentrations on the com- 
puted ozone concentrations for day 216 is shown in figure 3. All concentrations 
except those being varied are equal to the values listed in table II. In fig- 
ure 3(a), modeled values are shown for two cases, one (the upper curve) where 
the initial and upwind concentrations of formaldehyde and of the group "other 
aldehydes" were each equal to 10 ppb and the other (the lower modeled curve) 
where the concentrations of both species were equal to 0.5 ppb. The latter con- 
centrations were used to simulate all the days in this study. Figure 3(b) shows 
the effects of changing only the PIN 02 and RNO 2 concentrations from the initial 
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estimates of 10 ppb to the concentrations of 0*1 ppb which were used in this 
study and are given in table II. 

Measurements of concentrations of coupled species above the mixed layer 
Ci^top were not available; the concentrations were estimated (table II) and, 
with the exception of ozone, were held constant for all hours and days simu- 
lated- Estimates of other species concentrations above the mixed layer were 
obtained frcxn Duewer (ref. 14) when available and made independently when not 
available. Ozone received separate treatment from other species since higher 
ozone concentrations above the mixed layer may occur during air stagnation 
periods. When a subsidence inversion existed on a day being simulated, as was 
the case for day 205, the concentration of ozone above the mixed layer, the 
entrained concentration, was increased from its prescribed value of 0.025 ppm 
(ref- 14) to a higher value which was determined frcxn the magnitude of the mea- 
sured ozone curve for the previous day and from the diurnal variation in meteo- 
rological conditions. Previous studies (ref. 15) showed that the computed maxi- 
mum ozone concentration in this model is relatively insensitive to the entrained 
ozone concentration c^^top when ozone concentrations are measured at upwind 
stations. The rate of change of ozone concentration within the cell is influ- 
enced more by ozone advected into the cell than by vertical mixing. Upwind 
ozone advected into the cell, which is obtained from hourly measurements, 
includes a contribution frc«n entrained ozone. 

The RAPS data set included hourly averaged values for anthropogenic area 
and point-source emissions, Q in equation (1). Since emissions were defined 
in terms of a somewhat smaller grid than the cell size, the smaller grids were 
summed to obtain hourly emissions in the modeling region. Using measured mini- 
mum and maximum mixing-layer heights as input, appropriate plume-rise algorithms 
(refs. 16 and 17) were applied to hourly point-source emissions to account for 
emissions which might escape above the mixed layer. For the most part, the net 
effect of the plume-rise algorithms was to reduce emissions where low mixing 
heights were present. 

Emissions classes in the RAPS data set are nitrogen oxides (NOX) , total 
hydrocarbons (THC) , CO, nonreactive hydrocarbons, paraffins, total olefins, 
aromatics, and total aldehydes. Diurnal variations in the emissions of CO, NOX, 
and the reactive hydrocarbons (the sum of paraffins, olefins, arcanatics, and 
aldehydes) are shown in figure 4. Based on emissions data for Los Angeles 
(ref. 1 2) , emissions classes were partitioned into 90 percent NO and 1 0 percent 
NO 2 ; total aldehydes were divided into 66 percent formaldehyde and 34 percent 
other aldehydes; and total olefins were partitioned into 59 percent ethylene and 
41 percent other olefins. While the partitioning coefficients for St. Louis 
emissions may not have the same absolute percentages as those for Los Angeles, 
previous studies (ref. 15) show that small changes in emissions have little 
effect on simulated species concentrations. 

Using the RAPS emissions inventory (version available in Feb. 1 980) , emis- 
sions were compiled for days 205 (Friday) and 21 6 (Tuesday) . It was assumed 
that those weekday emissions were typical values for all summer weekdays since 
predicted oxidant levels were shown to be relatively insensitive to small vari- 
ations in CO and NO^ emissions (ref. 15). Therefore, emissions for day 216 were 
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also used to simulate day 226 (Friday ) , and emissions for day 205 were used to 
simulate days 1 59 (Monday) , 1 60 (Tuesday) , and 1 95 (Tuesday) . 


METEOROLOGY 

Meteorological data were measured at the surface at most RAM stations; 
the measurements used in the model were wind speed, wind direction, temperature, 
dew point, barometric pressure, and solar insolation. A single hourly averaged 
value for each parameter was formed using all stations within the cell 
(fig. 1). Wind speed and wind direction were used in two ways: to define trans- 

port through the cell using equation (1 ) , and to determine which stations upwind 
of the cell were used to specify concentrations ^i,bg advected into the cell. 
Dew-point measurements determined appropriate concentrations of water vapor. 
Temperature and barcmetric pressure measurements were used to maintain conserva- 
tion of mass in the model. Solar insolation was used to calculate UV-scaled 
photolytic rate constants. 

The height of the mixed layer characteristically increases during daylight 
hours (ref. 18), and the calculations are very sensitive to this growth param- 
eter. The rate of growth per hour varies with the time of day, with solar inso- 
lation, and with the larger-scale synoptic conditions. The minimum and maximum 
heights of the mixed layer vary with the day being simulated. Mixing heights 
were measured by the RAPS Upper Air Sounding Network (ref. 3) at only one 
station inside the cell, station 141, and at 6-hour intervals - approximately 
0400, 1000, 1600, and 2200 LST (table III). Measurements made at 0400 and 
1600 LST established the minimum and maximum heights of the mixed layer, respec- 
tively; measurements made at 1000 LST were used as guidelines for establishing 
the percentage growth for the rise of the mixed layer from the minimum value to 
a value at 1 030 LST. For the 0400 LST mixing-height measurement for 5 of the 
6 days in table III, the mixing depth as defined in the archived data tapes is 
listed as zero. This is due to the presence of a surface nocturnal inversion 
(ref. 19). In all five situations, however, a second low-level stable layer was 
present, and later in the morning it was eroded because of surface heating. In 
order to allow for a finite modeling volume at the start of a daily calculation, 
the top of the cell was set as the base of the elevated stable layer, even 
though it is not clear whether materials injected into the volume in the morning 
hours are mixed up to that base level. Knowledge of how quickly the low-level 
stable layer actually eroded would be helpful, but continuous measurements of 
the boundary-layer structure are needed for such an assessment. 

CO is used as an indicator of the growth rate of the mixed layer since 
changes in CO concentrations are derived primarily from volume changes within 
the cell and from transport through the cell rather than from chemical reactions 
(assuming that the emissions inventory adequately describes CO sources in the 
cell). When continuous measurements of the mixing-layer height are not avail- 
able, good agreement between measured and simulated CO concentrations can be 
obtained by adjusting the growth function for the rise of the mixed layer. This 
technique was used in model studies of the San Francisco Bay area using the 
Livermore Regional Air Quality (LIRAQ) model (ref. 20) . Using this approach, 
the early morning measured CO "traffic” peak was reproduced reasonably. Per- 
centage differences between the measured and modeled rates of growth are shown 
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in the last two columns of table 111 , and for two (195 and 205) of the 6 days 
(159f 160r 195, 205, 216, and 226), those differences were appreciable. The 
vertical temperature- gradient data near the ground indicated a strong surface 
inversion on the morning of day 205. For the remaining day (195), there was 
no early morning surface inversion, and the related meteorological parameters 
as well as the surface CO and NO air-quality data indicate good advection and 
surface heating in the morning. Surprisingly, however, the measured mixing 
height rose only 110 m (from 384 m to 494 m) from 0349 LST to 1048 LST, although 
it eventually grew to a level of 1905 m at 1550 LST. The mixing height of 
494 m at 1048 LST, if accurate, is therefore judged to have been ineffective in 
capping the volume, and this would explain the large discrepancy between the mea- 
sured and adjusted percentage growth rates for day 195. Therefore, when the 
discrepancies in the measured mixing heights for days 195 and 205 are considered, 
the estimated rates of growth of the mixed layer agree well with the measured 
growth rates. 

AS indicated in the preceding discussion, the variation of CO with time is 
used as an indicator of the growth of the cell volume due to the mixing-height 
rise of the transport and of the emissions components of the model. Only one 
chemical reaction involving CO, 


02 

CO + OH — ► CO 2 + HO 2 


occurs in the model, and it is of minor importance to the urban photochemistry. 
Since diffusion is neglected and only one-dimensiorial horizontal transport is 
considered in this model, the diurnal variation of CO concentration then becomes 
a good indicator of transport through the cell. The primary urban source of CO 
is vehicle emissions. If CO concentrations were dominated by vehicle emissions, 
then diurnal CO concentrations should closely follow fluctuations in traffic 
density. In particular, a 1- to 2-hour peak concentration should occur in the 
early morning followed by a gradual decrease in CO concentration (ref. 21). 
However, a sharp rather than a broad CO peak occurs in the measurements early 
in the morning, which reflects the beginning of the rise of the mixed layer, 
and thus the dilution of the CO concentration due to the expanding volume of the 
cell. 


In general, CO measurements are obtained at the surface close to the source 
emissions and, since CO concentrations are source dependent. They are spatially 
quite inhomogeneous. Likewise, vertical CO concentrations should be quite 
inhomogeneous and should decrease with altitude, especially during stable 
atmospheric situations. Simulations predict CO concentrations only for the mid- 
dle of the cell volume. Modeled CO concentrations should then be somewhat lower 
than measured concentrations. 

The starting times for the rapid growth of the mixed layer were determined 
by noting the times of the morning maxima in the CO curves. The CO peaks, and 
thus the starting times for rapid growth, occurred between 0630 and 0730 LST and 
depended upon the particular day being simulated. Table IV shows the hourly 
growth rate of the mixed layer used in the model. 


9 


L. 



The surface heating during the sunlight hours retains the base of the mixed 
layer at the surface. At sunset, the base of the mixed layer is no longer con- 
strained to the surface, and, in fact, may form a mixed layer aloft. Clean air 
then enters beneath the mixed layer, which produces discontinuities in surface 
concentration measurements. Since this physical process is absent in this model, 
agreement between measured and modeled species concentrations is expected to be 
poor after sunset. 


MODEL RESULTS AND DISCUSSION 

Air quality simulations were conducted for day 205 and day 21 6 to gain con- 
fidence in the use of the model. Simulations for these dates were shown previ- 
ously by Schere and Demerjian (ref. 4). 

Figures 5(a) to (d) show the averaged measured (circles) and modeled 
(squares) species concentrations within the cell for day 216 for CO, NO, NO 2 , 
and O 3 , respectively. The shaded areas represent a range of one standard devia- 
tion around the mean of species concentrations measured at all stations within 
the cell at each hour. Moderate wind speeds, which produce good advection 
through the cell, existed throughout the day. The simulations show good agree- 
ment with averaged measured concentrations. 

Changes in CO concentrations (fig. 5(a)) are derived primarily from volume 
changes in the cell, transport through the cell, and emissions inside the cell 
rather than frcan chemical reactions. For a given day, transport through the cell 
was determined by measurements of wind speed made at stations within the cell; 
it was assumed that CO sources within the cell were adequately described by the 
emissions inventory. Thus, the remaining variables that could be adjusted to 
provide better agreement between measured and simulated CO concentrations were 
the starting time for growth and the growth rate of the mixed layer, terms which 
describe volume changes in the cell. 

The sharp peak in the measured morning CO concentrations is an indicator 
of the expanding volume of the cell and, consequently, the starting time for the 
growth of the mixed layer. Thus, adjustments in the growth rate were used to 
obtain better agreement between measured and simulated CO concentrations. 

CO concentrations are source dependent with expected decreases in concen- 
tration as altitude increases. While CO measurements are made at the cell base, 
simulations predict concentrations in the middle of the cell volume; thus, pre- 
dicted CO concentrations should be less than measured concentrations. The growth 
rate of the mixed layer, particularly in the morning hours, was adjusted to 
reflect the decrease in predicted concentrations. The increase in concentra- 
tions of all species after 1 830 LST is expected. This increase is due to the 
breakup of the mixed layer at sunset, a process which is not described in this 
model. 

Diurnal variations in the concentration of NO, shown in figure 5(b), result 
from emissions, transport, chemical reactions, and volume changes in the cell. 
Transport occurs throughout the day. The concentration levels of NO in the 
early morning, from 0530 to 0630 LST, are determined by the magnitude of the 



emissions. Concentration levels of NO from 0630 to 1230 LST are dominated by 
emissions and increases in the cell volume; large fluctuations in the predicted 
concentrations of NO were obtained during this time period when the growth rate 
of the mixed layer and the emissions flux of NO were varied. After 1 230 LST, 
variations in the emissions flux and growth rate of the mixed layer had negli- 
gible effects on the concentration levels of NO, which indicated that chemical 
reactions dominated the levels of predicted species concentrations. Emissions 
of NO generally occur close to the ground, and measurements of species concen- 
trations are made at the ground; species concentrations are calculated for the 
middle of the cell volume. It is assumed that uniform mixing occurs in the 
model; however, when emissions are large, as they are in the early morning, 
instantaneous mixing of the emissions throughout the cell volume may not actu- 
ally occur. Therefore, when emissions are important contributors to the con- 
centrations of NO in the cell, from 0530 to 1 230 LST, predicted NO concentra- 
tions should be and are greater than measured concentrations. 

Two observations concerning the NO curves were made during the model vali- 
dation process and are worth noting. First, a more rapid decrease in NO concen 
trations from 0630 to 1 230 (obtained by increasing the growth rate of the mixed 
layer) increases the maximum ozone concentration predicted for the day. A prob 
able explanation for this observation is that NO scavenges ozone through the 
reaction 


O 3 + NO -► NO 2 + O 2 


Therefore, increased NO concentrations inhibit production of ozone. Secondly, 
the maximum ozone concentration occurs 4 to 5 hours after the crossing of the 
measured and modeled NO curves (at 1 230 LST) . The crossing (and thus the time 
of occurrence of the maximum ozone concentration) can be changed slightly by 
adjusting the growth rate of the mixed layer. 

The diurnal variations in the concentration of N 02 f shown in figure 5(c), 
result from emissions, transport, chemical reactions, and volume changes in the 
cell. It is not clear frcwi examining the fluctuations in the averages of mea- 
surements made inside the cell whether a single process dominates at any given 
time. The shape of the simulated curve and the magnitudes of simulated NO 2 con 
centrations should agree with those obtained from measurements in order to 
obtain good agreement between measured and modeled NO and O 3 curves. 

Changes in O 3 concentrations for day 216 (fig. 5(d)) are attributed to 
transport and chemical reactions. Ozone is scavenged at the ground by surface 
reactions and dry deposition, both of which are neglected in this model. As a 
direct consequence, modeled ozone concentrations for afternoon hours, which 
represent concentrations in the middle of the cell, are expected to be greater 
than concentrations measured at the surface. A disparity in ozone concentra- 
tions of the order of 1 5 percent was observed by Duewer et al. (ref. 20) in 
their measured and modeled (LIRAQ) results for the San Francisco Bay area 
stations. The simulated O 3 concentrations obtained in this work show good 
agreement with average measured concentrations. The modeled ozone maximum 
of 121 .5 ppb is 11.7 percent higher than the measured value of 108.7 ppb; 
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both the position of the maximum and the shape of the simulated curve agree well 
with the measured curve. 

Predicted ozone concentrations in the morning hours (0530 to 1230 LST) are 
consistently lower than measured concentrations. Although it is not clear why 
the predicted concentrations are underestimated, one possible explanation is that 
the chemical mechanism overcompensates for ozone scavenging by NO. This dis- 
crepancy also appears in simulations of other days and certainly warrants study 
in the future. 

An inversion existed on day 205 (simulations shown in fig. 6) which pro- 
vided quite different meteorology from day 216. Ozone from the previous day was 
trapped aloft since the wind speeds were insufficient to introduce clean air 
above the urban area. Consequently, an elevated entrainment concentration of 
60 ppb ozone was used (ref. 4). Since previous studies (ref. 15) showed that 
the model results are relatively insensitive to the concentration of ozone aloft 
when measured concentrations at stations upwind of the cell are available, the 
entrainment concentration of 60 ppb for day 205 compared with the "clean air" 
value of 25 ppb indicates only the existence of an inversion but not necessarily 
the exact strength of the inversion. The measured morning growth rate for the 
rise of the mixed layer was less than the growth rate for day 216. The esti- 
mated morning growth rate for day 205 was similarly reduced relative to day 216 
in order to decrease the computed maximum ozone concentration and increase the 
width of the ccxnputed ozone-concentration curve. The computed ozone maximum had 
a value 13.7 percent higher (94.5 ppb) than the measured ozone concentration 
(83.1 ppb). 

Simulations for days 205 and 216 were performed to parameterize the par- 
titioning coefficients for total hydrocarbon concentrations into hydrocarbon 
species and to determine appropriate initial concentrations of other species; 
verifications of the parameter izations were obtained through additional simu- 
lations of days 159, 160, 195, and 226. A critical variable in the simulations 
was the growth rate of the mixed layer. When days 159, 160, 195, and 226 were 
simulated, modeled CO concentrations were larger than measured concentrations 
throughout the day, which suggests that representative CO emissions were too 
large in the early morning. The introduction of a small early morning growth 
rate (less than 2.5 percent per hour) prior to the morning peak in the measured 
CO curve allowed modeled CO concentrations to remain lower than measured results 
throughout the day. 

Simulations for day 160 are shown in figure 7. The vertical temperature 
gradient measured between the surface and the top of a 30-m tower was negative, 
which indicated that no inversion was present. Ozone reached a maximum measured 
concentration (averaged over all stations inside the cell) of 124.5 ppb on the 
previous day. The measured concentrations of CO, NO, and NO 2 in the morning 
hours were much higher than the morning concentrations of those species on previ- 
ous days (days 205 and 216); in fact, the measured average concentration of NO 
at 0530 LST, the starting time for the simulation, was about three times as 
large as initial NO concentrations for days 205 and 216. The increased initial 
concentrations for this day could introduce difficulties in the simulation, 
since the model was validated using the lower concentrations of days 205 and 216. 
A distinct morning peak in the measured CO and NO concentrations (which indi- 



cates the starting time for the growth of the mixed layer) was not evident. 
Although fairly large increases in the concentration for this day were present/ 
no difficulties were encountered in simulating concentrations. Good agreement 
with the measured concentrations was obtained in the simulations of CO, NO, 
and NO 2 . The maximum in modeled ozone concentrations of 174 ppb also compared 
quite favorably with the measured maximum concentration of 1 63 ppb (a difference 
of 7 percent) . 

Figure 8 shows simulations for day 195. Moderate wind speeds measured for 
for this day provided good advection through the cell. Initial concentrations 
of 00, NO, and NO 2 were reduced relative to those for the validation days, but 
the initial concentration of ozone was at least twice as large. These variations 
in starting concentrations relative to the validation days provided a further 
test of the applicability of the model to wide ranges of concentrations. While 
the modeled curve for CO showed particularly good agreement with the measured 
CO curve in both shape and magnitude, generally good agreement was obtained 
between all measured and modeled curves. There was no difference in magnitude 
between the measured and modeled ozone-concentration maxima. 

Figure 9 shows the simulations for day 159. The wind speed was very low, 
and the measured diurnal ozone curve was broad and flat, which indicated high 
air pollution potential. The modeled ozone concentration reached a maximum 
value of 142.5 ppb at 1530 LST, which was 16.6 percent higher than the measured 
value of 124.5 ppb which occurred at 1330 LST. Good agreement was found between 
the modeled and measured concentrations of CO and NO; poorer agreement was 
obtained for the NO 2 curves in the morning hours. The discrepancies obtained 
for the NO 2 curves compared with the CO and NO curves are not surprising since 
a major source of NO 2 is industrial emissions but the major source of CO and NO 
emissions is vehicular activity. Thus, although the representative emissions 
used in the model may adequately represent vehicular emissions, variations in 
industrial emissions may account for the variations in the NO 2 curves for this 
day. 


Maximum ozone concentrations for day 226 (fig. 10) were quite high. Cloudy 
sky conditions existed on the previous day with a relatively small buildup of 
ozone. The trailing edge of a high-pressure system was in close proximity to 
the simulation area. While the measured ozone curve was quite broad, adjust- 
ments in the growth rate of the mixing-layer height to produce reasonable 
agreement for CO, NO, and NO 2 simulations produced a narrow simulated ozone 
curve with a maximum concentration that was 3 percent larger than the measured 
peak concentration. 

Table V summarizes the model performance for predicting the time of occur- 
rence and magnitude of the ozone-concentration maxima. Modeled ozone maxima 
are greater than measured maxima by an average of 8.7 percent, a reasonable 
increase and well within the 1 5 percent maximum increase noted by Duewer et al. 
(ref. 20) in studies of the San Francisco Bay area using a more complex (LIRAQ) 
model. The time of occurrence of the maximum concentration of ozone was not 
simulated as accurately, possibly a consequence of the simplified treatment of 
chemistry in the model. Likewise, for all data modeled here, the predicted con- 
centrations of ozone in the early morning hours were consistently underestimated, 
perhaps a consequence of over compensating for the morning scavenging of O 3 by 



NO in the chemical mechanism. The good agreement in predicting the ozone maxima 
and the diurnal variations in other species concentrations can be used as one 
criterion to assess the usefulness of the model. 

The use of the model should be undertaken with some caution since the pro- 
cesses of dry deposition, surface scavenging, and heterogeneous reactions have 
not been included. Additionally, although it is assumed in the model that 
species concentrations are homogeneous with respect to altitude, variations in 
concentrations with respect to altitude have been measured for other urban areas 
(ref. 22). However, the model performs quite well with data from the limited 
number of parameters usually measured in an air-quality experiment. The addi- 
tion of more detailed physics into the model to describe the urban troposphere 
is warranted only when data describing additional physical processes are 
available. 


SUMMARY OF RESULTS AND CONCLUDING REMARKS 

A single-cell Eulerian urban-scale air-quality simulation model was devel- 
oped and validated using six clear-sky data sets obtained during the 1 976 EPA 
Regional Air Pollution Study (RAPS) in St. Louis. The urban area is described 
in the model by a single cell with a fixed horizontal length and a variable ver- 
tical height which depends on the height of the mixed layer. Other parameters 
input to the model were wind speeds, an emissions inventory, chemical mechanism, 
and measured species concentrations upwind, inside, and above the cell. 

Predicted species concentrations were relatively insensitive to small vari- 
ations in CO and NO^ emissions and to the concentrations of species above the 
cell top which are entrained as the mixed layer rises. A greater sensitivity 
to other input parameters was shown. In particular, nonzero concentrations of 
species inside and upwind of the cell were required to initiate the chemical 
mechanism. Because of the rather small horizontal dimensions of the single 
cell (20 km by 20 km) , predicted species concentrations were sensitive to the 
advection of upwind ozone and its precursors. Chemical reactions were domi- 
nated by the chemistry of the reactive hydrocarbons. Two key variables in the 
model had to be estimated, the coefficients used to partition emissions into 
species concentrations and the hourly growth rates of the mixed layer. Since 
chemical reaction rates are concentration dependent, simulated ozone concentra- 
tions were highly dependent upon the magnitude of the scale factor used to par- 
tition the concentrations of reactive hydrocarbons into hydrocarbon species. 

For 4 of the 6 days simulated in this study, the estimated growth rates of .the 
mixed layer agreed quite well with measured values. Significant changes in the 
predicted concentrations of ozone were also obtained when photolytic reaction 
rates were scaled from theoretical to atmospheric values to reflect the attenua- 
tion of solar radiation in the atmosphere by particulates and aerosols. Infor- 
mation on the actual concentrations of ozone precursors and on the homogeneity 
in species concentrations and meteorological parameters in the modeling volume 
is also needed to accurately predict ozone concentrations in the middle of the 
cell volume. 

Good agreement was obtained between measured and modeled concentrations of 
NO, CO, and NO 2 for all days simulated. The concentrations of ozone in the 



early morning hours were consistently underestimated, a probable consequence of 
overccxnpensating for the scavenging of ozone by NO in the chemical mechanism; 
however, good results were obtained in predicting the magnitudes of the ozone 
maxima. The predicted maximum concentrations of ozone for the six simulations 
were an average of 8.7 percent higher than the measured maximum concentration. 
This is an anticipated and reasonable difference considering that the model 
does not treat surface scavenging and dry deposition, processes which could 
decrease the ozone concentrations that were measured at the ground. The most 
difficult data sets to simulate were those where a strong inversion was present; 
simulations of these data sets also showed good agreement with the measured data 
sets. 


The model formulation could be improved if additional types of data were 
measured. Sufficient chemical- kinetic data are not available to include aro- 
matic and heterogeneous chemical reactions in the chemical mechanism. Surface 
scavenging and dry deposition deplete measured species concentrations at the 
surface; measurements of the types and quantities of particles lost from dry 
deposition are needed. Species concentration profiles would allow for both the 
parameterization of the effect of surface scavenging and for the assessment of 
how surface scavenging and dry deposition contribute to changes in species con- 
centrations inside the cell. Surface scavenging and dry deposition may have 
negligible contributions to the changes in species concentrations in the center 
of the cell and thus to the predictive capability of the model in the center of 
the cell; however, the contributions of these processes cannot be judged using 
the surface-station data in the RAPS data set. The parameterization of the 
cessation of surface heating and the subsequent formation of a mixed layer aloft 
would require much more information about ccmiplex meteorological processes and 
species concentrations. Even given the improved data set, the chemical mech- 
anism (which was set up for daytime chemistry) would require considerable test- 
ing and perhaps modification. In any event, the parameterization is beyond the 
scope of the purpose of this model. 

The model performance and validation could also be enhanced if several addi- 
tional types of data were measured. Periodic measurements of concentrations of 
reactive hydrocarbon species for validation purposes and hourly measurements of 
the height of the mixed layer would eliminate the present need to estimate these 
quantities. Since transport through the center of the cell is a direct function 
of wind speed, the addition of vertical wind profiles rather than the use of sur- 
face winds would permit a more accurate treatment of transport through the cell. 
The addition of vertical concentration profiles of the criteria pollutants would 
give a better estimate of how well the model performs, since species concentra- 
tions are calculated in the center of the cell volume. Data in some of these 
areas were obtained during the RAPS field program, but none were archived for 
general use during the time of this work. 

With the possible requirement of adding local anomalies, the model should 
be quite easily adapted for use in other urban centers. Measurements of species 
concentrations of reactive hydrocarbons are needed, and the emissions inventory, 
which is unique for each urban center, should be changed. Adjustments in the 
parameter izations for the growth of the mixed layer and for the scaling factor 
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for the photolysis of NO 2 are necessary. Initial species concentrations and 
the concentrations of background and entrained species might also require 
modification. 

The single-cell model can be used with success when sufficient groundbased 
measurements are available. Since the surface-station measurements used to 
validate this model are of the type normally available in any large urban area 
from Air Pollution Control Board monitoring stations, this type of model is 
more readily adapted to other urban areas than more canplex models requiring 
more detailed data sets. Still, this model requires significant amounts of 
data to initialize and simulate air quality for a given day. Unless accurate 
and periodic measurements of the height of the mixed layer are available, the 
model cannot fully characterize pollution episodes. 

Given the limitations in the data set used for this model validation 
study, the present single-cell photochemical air quality model is a suitable 
tool for studying the chemistry and transport in an urban atmosphere. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 11, 1 981 
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TABLE I 


Hydrocarbon class 


Nonreactive 

hydrocarbons 

Alkanes 

Ethylene 

Other alkenes 


Aromatics 


AVERAGE HYDROCARBON CONCENTRATION LEVELS 
FOR LOS ANGELES AND ST. LOUIS 


Los Angeles (1 965) 

St. Louis 
(present study) 

Magnitude, 

ppm 

THC, percent 

Scale 

THC, percent 

3.22 

84.99 

1 .1 31 2 

96.1 45 

.341 2 

9.01 

.2830 

2.55 

.06 

1 .58 ’ 

.2278 

.36 

.042 

1 .11 

.2838 

.31 5 

.085 

2.24 

.281 2 

.63 


9 




TABLE II.- CONCENTRATIONS OF COUPLED SPECIES 


Species 

Concentration, ppbv* 

^initial 

^bg 

^top 

C2H4 

0.4 

0.4 

0.4 

Other olefins 

.4 

.4 

.4 

H2C0 

.5 

.5 

.5 

Other aldehydes 

.5 

.5 

.5 

03 

+m 

m 

^25.0 

NO 

m 

m 

1 0.0 

NO2 

m 

m 

2.0 

CO 

m 

m 

1 20.0 

Alkanes 

35.0 

35.0 

35.0 

HNO2 

.1 

.1 

.1 

RNO2 

.1 

.1 

.1 

HNO3 

1 .0 

1 .0 

1.0 

RNO3 

.001 

.0 

.0 

H2O2 

.01 

.01 

.01 

PAN§ 

.001 

.0 

.0 


*°initial “ Initial concentration inside the cell at 0530 LST; 
cjjg = Background concentration, concentration upwind of cell; 
c^op = Entrained concentration, concentration above mixed layer. 
"^Measured hourly. 

$For day 205, <^ 3 , top =60.0 ppbv. 

^Peroxyacetyl nitrate. 








TABLE IV.- HOURLY GROWTH RATE OF 


MIXED LAYER USED IN MODEL 


Day 

Time, 

LST 

Percentage 

growth 

1 59 

0530-0630 

1 .5 


0630-0730 

7.0 


0730-1 030 

47.1 


1 030-1 630 

44.4 

1 60 

0530-0630 

2.0 


0630-0730 

6.0 


0730-0830 

40.0 


0830-1 030 

32.5 


1 030-1630 

19.5 

1 95 

0630-0730 

2.0 


0730-0830 

40.0 


0830-0930 

20.0 


0930-1 530 

38.0 

205 

0630-0730 

5.0 


0730-1 030 

23.0 


1 030-1 630 

72.0 

21 6 

0630-1 030 

64.0 


1 030-1 230 

16.0 


1 230-1 530 

20.0 

226 

0530-0730 

4.5 


0730-0830 

50.0 


0830-0930 

18.0 


0930-1 030 

1 5.5 


1 030-1 630 

1 2.0 

^ 
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